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Abstract. Recent intensive observations of the BL Lac object OJ 287 raise a lot of questions on the models 
of binary black holes, processing jets, rotating helical jets and thermal inst abilit y of slim accretion disks. After 
carefully analyzing their radio flux and polarization data, Valtaoja et al. (2000) propose a new binary model. 
Based on the black hole mass of 4 x 10 s Mq estimated with the tight correlations of the black hole masses and the 
bulge luminosity or central velocity dispersion of host galaxies, we computed the physical parameters of the new 
binary scenario. The impact of the secondary on the accretion disk around the primary black hole causes strong 
shocks propagating inwards and outwards, whose arrival at the jet roots is identified with the rapid increase of 
optical polarization and the large change of polarization angle at about 0.30 yr after the first main optical flare. 
An increase of optical polarization, a large rotation of positional angle and a small synchrotron flare at 2007.05 
between the optical outbursts at 2006.75 and 2007.89 are expected by the model. With the estimated parameters, 
we predicated an increase of 7-ray flux appearing about 5 days after the first optical/IR peak, which is consistent 
with the EGRET observations. 
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1. Introduction 

OJ 287 (0851+202) at red-shift z = 0.306 is one of the 
most active and best studied B L Lac objects. While it was 
identified in 1967 (Dickel et al. 1967 ), its optical observa- 
tions were dated back to 1890s, which makes it be one of 
the best candidates for searching for long-term periodic- 
ity. Ba sed o n the historical optical light curves, Sillanpaa 
et al. (|1988|) found the large optical outbursts of OJ287 
recurrent with a period of 11.65 yr and proposed a binary 
black hole model, which predicted an outburst in the 1994 
fall. 

Although the predicted large optical outbursts were 
detected in 1994 and 1995, it was found that the out- 
bursts are double-peaked and the binary black hole model 
has to be modified (Sillanpaa et al. 1996 ). Among the pro- 
posed models to explain the periodicity of O J 287, those 
of binary black holes (Lehto & Valtonen 1996), procession 
jet (Katz |l997p , rotating helical jets (Villata et al. |998|) 
and the thermal instability of slim accretion disk (Liu et 
al. 1995) can explain the periodicity in optical band well. 
However, radio and polarization observations show that 
the first of the two peaks is narrow and thermal, and 



the second one is broad and synchrotronic (Valtaoja et 
al. 2000; Pursimo et al. 200C). Both flares should be syn- 



chrotronic and identical in the processing/helical jet mod- 
els while thermal and similar to each other in the binary 



black hole model (Valtaoja et al. 200C). Another difficulty 
with the binary model is that it requires an extremely 
massive primary of M = 1.5 x 10 10 M o (Pietila |1999| ), 
which is at least one order of magnitude larger than the 
super-massiv e blac k hole ma sses o f other BL La c obje cts 
(Barth et al. |2002| ; Wu et al. [2002| ; Falomo et al. |2002|) . 

In the binary black hole model, the outburst ther- 
mal optical/IR photons should be up-scattered into 7-ray 
regime by the relativistic electrons in the jets via inverse 
Compton scattering process. From the parameters given 
by Pietila (1999), the estimated time for photons to travel 
from the impact site to the jets is about 20 days and im- 
plies a delay of about three weeks between the optical/IR 
outbursts and high energy 7-ray radiations. Such a large 
lag is not consistent with the EGRET 7-ray observations 



(Shrader et al. |1996j). 

To resolve the difficulties of the models mentioned 
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above, Valtaoja et al. (2000) suggest a new binary scenario 
to explain the observations of radio flux and polarization. 
In their model, the first thermal flare of the outbursts is 
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caused when a secondary black hole hyper-sonically passes 
through the accretion disk around the primary, and the 
second synchrotron flare after a viscous time scale i V j S is 
from the jet as the perigee passage enhances the accre- 
tion rate in the accretion disk, leading to increased jet 
mass flux and formation of strong shocks down the jet. As 
for such a scenario a constant period agrees with all the 
sharp outbursts (the first peak), a precession of the binary 
orbit and an extreme primary black hole are not needed 
(Valtaoja et al. |2000p . However, as the lack of accurate 
knowledge of the primary black hole mass, Valtaoja et al. 
( 2000D did not estimate the parameters of the new system. 
In this letter, we extend Valtaoja et al's model to explain 
all the observations and try to estimate the physical pa- 
rameters of the new scenario. 

After a brief summary of the observations in Sec. [| we 
describe the model in Sec. |. In Sec. |, we try to estimate 
the parameters of the model. Our discussions and conclu- 
sions are presented in Sec. [|. Hq — 50kms _1 Mpc _1 and 
qo = 0.5 are used throughout the paper. 



2. Observations 

The optical outbu rsts o f OJ287 have a period of P = 
11.886 yr (Kidger [2000D and are do uble- peaked with a 
spacing of 416 days (Sillanpaa et al. |1996[ Kidger |2000D . 
The first flare is thermal and narrow without correspond- 
ing radio outburst and the second one is non-thermal and 
broad with strong outburst in radio (Valtaoja et al. 2000; 
Pursimo et al. 2000 and references therein). The variabil- 
ity in the IR band is similar to that in the optical band 
(Kidger et al. |1995| ; Pursimo et al. |2000D . The first flare 
in the nearest 1994-1995 outbursts peaked at JD1994.86 
and the second one did at 1996.00. A small UV flare was 
detected during the first optical flare and no observation 
was made around the second one (Pursimo et al. 200C). 

The observations made around JD1994.88 showed no 
evidence for any X-ray outburst, but weak flux variations 
of time-scale as short as the time limit of the instrument 



(96 m inute s) were detected (Idesawa et al. 1997 Pursimo 
et al. 2000| ). EGRET observations integrating from 10 to 
15 November 1994 showed that in the 7-ray band the ob- 
ject was three times brighter than earlier, when it was in a 
low optical state (Shrader et al. 1996; Pursimo et al. 2000) 

The optical polarization is very low during the first 
flare but very high during the second flare. A substan- 
tial increase of optical polarization and a large change of 
positional angle within one day at JD199 5.16 have been 
obser ved (Smith et al. 1987 ; Takalo 1994 ; Pursimo et al. 
2000| ; Efimov et al. |2002| ). The variability of radio polar- 
ization is similar but with a time delay to that of optical 
polarization (Valtaoja et al. 2000 ; Pursimo et al. 2000 ). 

The host galaxy of OJ 287 has been marginally re- 
solved recently and has an absolute magnitude of A/r = 



3. A revised binary black hole model 

The new binary scenario initially proposed by Valtaoja et 
al. ( 2000 ) and described in the introduction can be ex- 
tended to explain the optical polarization, the X-ray vari- 
ation and the 7-ray radiation of OJ287. 

In addition to an enhancement of the accretion rate, 
the perigee passage also causes strong shocks in the accre- 
tion disk, which propagate inwards and outwards at the 
speed of sound. The shocks arrive at the jet roots after 
a sound traveling time scale At s and cause cause a very 
small increase of the jet flow and shocks down the jets, 
compressing the magnetic fields in the jets and therefore 
causing a rapid increase of the optical polarization and 
large rotation of positional angle about At s ~ 0.30 yr af- 
ter the first optical flare. 

The thermal optical/IR photons of the first main flare 
feed the jets of the primary black hole and are up-scattered 
into the 7-ray regime by the relativistic electrons in the 
jets via inverse Compton scattering. The 7-ray outbursts, 
which have a delay of light traveling time from the impact 
site to the jets, may have been detected by EGRET. The 
detected weak X-ray variability may be generated around 
the secondary black hole, as the secondary is also an X-ray 
radiation source. 

The apogee passage of the secondary black hole 
through the accretion disk may be too weak and unde- 
tectable. 



4. Parameters of the model 

4.1. Masses of the black holes 

It was found recently that the masses of the center black 
holes in active and inactive galaxies follow tight correla- 
tions with the bulge velocity dis persion s (e.g. Gebhardt 



et al. |2000| ; Merritt & Ferrarese [2001a[ |2001b| ) and the 
bulge lu mino sities (e.g. Magorrian et al. 1998j ; McLure & 
Dunlop 2002 ). The host galaxy of OJ287 was marginally 
resolved only recently of an effective radius r c = 0.72" 
and an absolute magnitude Mr = —23.23 (Heidt et al. 
1999Q . A black hole mass of M ~ 4.6 x 1O 8 M is obtained 
with the absolute magnitude together with the correlation 
of black hole mass and bulge luminosity (M — Ltmigc) for 
active galaxies given by McLure & Dunlop ( 2002 ). With 
t he m orphological data of the host galaxy in Heidt et al. 
( |1999| ), we estimated the central velocity dispersion er, us- 
ing the fundamental plane of a tight correlation of velocity 
dispersion, effective radius and average surface brightness 
for elliptical galaxies (Bettoni et al. 2001). A second value 
of the black hole mass M ~ 3.2 x 1(PM^ was derived with 
the estimated velocity dispersion er and the M — a corre- 



-23.23 (Heidt et al. 1999) 



lation given by Merritt & Ferrarese (20011). The two esti- 
mated values of the black hole mass agree well with each 
other and are consistent with the upper limit M ^ 10 9 M Q 
obtained by Valtaoja et al. (2000). We used their mean 
value M = 4 x 10 8 M Q in the following calculation. 
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For the secondary, one can set only an upper limit 
m/M® ^ 10 8 , if the variability of the X-ray flux is from 
the accreting gas around it. 

4.2. Standard thin accretion disk or ADAF? 

While it has been suggested that the ac cretio n disk in BL 
Lac ob jects is optically thin (Rees et al. 1982 ; Cavaliere & 
D'Elia pOOl ), it is still possible that OJ 287 is exceptional 
and has an optically thick and geometrically thin accretion 
disk (e.g. Sundelius et al. |l997j Pietila et al. f999h. 



For an accretion disk, the viscous time-scale is t V i S ~ 
R/v ~ R 2 /i> — Rt s /(aH) where t s — R/c s . If we take 
i V is — 1-14/(1 + z)yr and t s ~ 0.30/(1 + z) yr, we have 
H/R <; 0.3 for a ^ 1, which implies a geometrically thick 
disk. If we assume that large rotation of optical polariza- 
tion angle (PA) around the time of the first optical out- 
burst is caused by the passage of the secondary through 
the accretion disk, the nearly simultaneity between the PA 
rotation and optical outburst implies that the outburst ra- 
diation region becomes optically thin immediately when it 
is generated by the passage. So, it is naturally to assume 
the accretion disk being optically thin and geometrically 
thick, i.e. an advection-dominated accretion flow (ADAF). 

For an ADAF, the self-similar solutions of sou nd speed 
c s and radial velocity v are (Narayan et al. 2000 ), respec- 
tively, c s = co«k and v — vqVk, where vk is the Keplerian 
velocity, cq = (2/5) 1 / 2 and vo — 3a/ '5 for the politropic 
index 7 = 5/3. Therefore, we have 



10 L 



3 U 



0.28, 



(1) 



which is very close to the favored value a 
erature (Narayan et al 



1998) 



0.25 in lit- 
The value of a slightly de- 
pends on the polytropic index 7 and we have a ~ 0.31 for 
7 = 7/5. 

4.3. Parameters of the binary system 

If the secondary black hole is comparable to the primary, 
the primary would rotate around their mass center with 
a large amplitude, which leads up to in a significant si- 
nusoidal twist of jets. Such a twist was not detected by 
VLBI (e.g. Tateyama et al. |l996t Roberts et al. [19871) . So, 
we have m <C M and the major axis a of the secondary 
orbit is 



1 

a — - 



Pc 



2 \-ktg (1 + z) 



2/3 



r G = 4.1 x 10 2 A/ 8 2/3 r G 



(2) 



where r G is the Schwarzschild radius of the primary and 
M s = Mj (4 x 10 8 M Q ). 

The pericenter passage of the orbit is at radius 



R„ 



c s t s /(l +z) = 88M 8 2/3 r G . 



(3) 



If the accretion disk is a standard thin high-a disk, the 
orbit of the secondary may be coplanar with the accretion 
disk due to disk-secondary black hole interaction (Ivanov 



et al. 1995 ; cf. Valtaoja et al. 2000 ). However, the situa- 
tion is unclear for an advection-dominated accretion flow, 
as the interaction between the secondary and the disk is 
negligible (Narayan 2000) and the coplanar time scale is 
of order of viscous time scale or larger (Ivanov et al. 19! 
Scheuer & Feiler 1996). In the following discussion, we as- 



sumed that the orbit of the secondary is roughly coplanar 
with the accretion disk and the secondary passes through 
the accretion disc twice per orbit. Then, the eccentricity 
of the orbit is e ~ 0.78 and the apogee is at radius 



Pu 



7.2 x 10 2 M 8 2/3 



r G . 



(4) 



The IR/optical photons, which are up-scattered into 7- 
ray regime by the relativistic electrons in the jets, travel 
from the pericenter passage to the jets within a time scale 
Aiy ~ Rmin/c = 4.0M 8 days. So, the model predicts an 
increase of 7-ray flux with a lag of (1 + z)A£ 7 ~ 5 days 
after the main flare, which i s cons istent with the EGR ET 
observations (Shrader et al. 1996; Pursimo et al. 2000] ). 



4.4. The apogee passage 



From the self-similar solutions (Narayan et al. 2000 ), the 
ratio of the mass densities between the two passages is 
Pmax/Pmin ^ (-Rmin/-Rma X ) 3/2 = 4.3 x 10~ 2 . As the radia- 
tion of optically thin hot gas is bremsstrahlung and pro- 
portional to the density square, the radiation is at least 
three order of magnitudes fainter at the apogee than at the 
pericenter. The difference of temperature, T max /T m j n ~ 
(-Rmin/-Rmax) = 1-2 x 10 -1 , makes the radiation at i? max 
even weaker and in lower frequency. Therefore, the burst 
at the apogee is undetectable. 



5. Discussions and Conclusions 

Basing on the observations of O J 287 in multi-frequencies 
and of the host galaxy, we extended and estimated the pa- 
rameters of the revised binary model initially suggested by 
Valtaoja et al. ( |2000|) . In the model, the mass of the pri- 
mary black hole is 4 x 10 8 Mq, which is consistent with the 



upper limit obtained by Valtaoja et al. (2000) and about 
40 times smaller than what the L&V's model has given 
(Pietila 1999). As the semi-major axis is, respectively, 
about 30r G in the L&V's model and about 4.1 x 10 2 r G in 
the revised model, the two black holes are about ten times 
closer in the L&V's model than in ours (in units of the 
Schwarzschild radius). With such a less massive primary 
black hole and larger semi-major axis, the gravitational 
radiation is much weaker in the new scenario, though we 
cannot exactly estimate its value as we have only an upper 
limit of the companion mass. It can be expected that the 
lifetime of the binary system is much longer in the new 
model than that in the L&V's which is only about 10 4 yr. 

As we took the constant period of P — 11.886 yr as in 
Kidger ( [2000 ), the predicted next outbursts (Kidger 2000) 
are that the first thermal optical flare is at 2006.75 (1 Oct 
2006) and the second synchrotron flare is at about 2007.89 
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(21 Nov 2007). The predictions are consistent within the 
uncertainties of the period with those given by Valtaoja 
et al. ( 2000| ). In addition, an increase of optical polariza- 
tion and large rotation of PA at 2007.05 (18 Jan, 2007) 
and brightening in 7-ray on 5 Oct, 2006 (at 2006.76) are 
expected by our model. 

We identify the small X-ray variations observed about 
one week after the first flare with radiations from the hot 
gas around the secondary black hole. However, it is possi- 
ble that the X-ray variability is not relevant to the peri- 
center passage event and cannot be detected during the 
2006-2007 outbursts. Another possibility is that the arrival 
signal of sound wave traveling from the impact site to the 
jets is the variation in X-ray continuum instead of optical 
polarization. If so, an extremely high eccentricity e = 0.98 
is needed. The orbit of the secondary in such a system is 
completely embedded in the thick accretion disk and no 
outburst can be observed. Any binary system with such 
an extremely-high initial eccentricity would very rapidly 
reduce its eccentricity to e 0.8 due to gravitational radi- 
ation and the interaction between the secondary and the 
accretion disk (cf Artymowicz 1992). 

If the traveling time scale of shock waves is as short as 
i s /2, the eccentricity of the binary system is e ~ 0.87 and 
most part of the orbit must be embedded in the accre- 
tion disk. No sharp optical outbursts can be detected. If 
the traveling time scale of the shock waves is longer than 
t s , e.g. 2i s , the orbit becomes very circular with e = 0.66 
and the two passages are very close to each other and 
Rmin/ Rmax — 0.21. The outbursts caused by the apogee 
passage become significant and detectable. In this case, an 
intersection of two periodic variabilities of p = 11.886 yr 
should have been observed. One possibility is that the or- 
bit is circular of a ~ 640rQ (Note that the observed period 
fobs = -Porbit/2) and we have actually observed two equal 
passages with i? m i n ~ R m ax — a- In this case, the traveling 
time scale of shock waves is about 7.3i s ~ 2.2 yr, which 
implies an unreasonably large a ~ 2.0. Our conclusion is 
that the shock traveling time scale of shock waves is be- 
tween about i s /2 and 2i s . It is reasonable to assume that 
the rapid increase of the optical polarization and the large 
rotation of positional angle after about 0.30 yr following 
the first flare are caused by the shock waves arriving at the 
jet foot-points. As the shock waves may increase the mass 
flow into jet and cause shocks down the jets, a small syn- 
chrotron outburst corresponding to the change of optical 
polarization is also expected. 

We have assumed that the orbit of the secondary is 
roughly coplanar with the accretion disk in the calculation 
of the eccentricity e and i? max - If they are not roughly 
coplanar with each other, the system would have a much 
higher eccentricity and closer apogee. This problem has 
already been discussed before. 

We used the self-similar solutions to a station- 
ary and vertically averaged one-dimension ADAF. Two- 
dimensional numerical simulation (Igumenshchev & 



Abramowicz 2000) showed that an ADAF with higher a 



(<; 1) is convectively unstable. For an ADAF with lower 



viscosity (a ^ 0.1), no bipolar or large scale circular mo- 
tion is present and i v i s should be passage- independent. 
For ADAF with a moderate viscosity a ~ 0.3, while 
a self-similar solution is adequate, time-dependent bipo- 
lar outflows and large scale circulations are important 
(Igumenshchev & Abramowicz 2000| ). Thus, a passage- 
dependent time scale i v ; s is not inconsistent with our 
model. However, to fully understand the time dependences 
of the inner accretion disk, two or even three dimensional 
numerical simulations are needed. In our model, a time- 
dependent ivis gives a time-dependent viscous parameter 
a. Other estimated parameters are independent of i v i s . 
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